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Abstract. We have developed a dynamical model for a unified description of the pion-nucleus scattering and photo- and 
neutrino-induced coherent pion production on nuclei. Our approach is based on a combined use of the Sato-Lee model for the 
electroweak pion production on a single nucleon and the A-hole model of pion-nucleus scattering. Numerical calculations are 
carried out for the case of the '^C target. After testing our model with the use of the pion photo-production data, we confront 
our predictions of the neutrino-induced coherent pion production reactions with the recent data from K2K and MiniBooNE. 
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INTRODUCTION 

In recent years both the theoretical and experimental studies on the coherent pion production in neutrino-nucleus 
reactions in few-GeV region have been intense. These studies were originally inspked by the on-going and future 
neutrino experiments. Some interesting results from recent experiments have further encouraged the theoretical studies 
of this process. Namely, the K2Kfl'] and SciBooNEfJ] collaborations reported no evidence of the charged-current (CC) 
process, while the MiniBooNE collaboration found a signal of the neutral-current (NC) process. [3] This discrepancy 
must be resolved since the isospin consideration predicts the simple relation Occ ^ ^(^nc- The on-going data analysis 
of the coherent process have been trying to provide more detailed information including the muon and pion kinematics 
for the CC process, and some preliminary results have been presented in this workshopl4, 5J. Considering these 
experimental developments, it is important to develop a rigorous theoretical model for analyzing the available and the 
forthcoming data. This is the objective of our work. 

There are mainly two types of theoretical approaches to the coherent pion production, i.e., a PCAC-based model 
and a microscopic model. We pursue the latter option. In a microscopic approach, the starting point is a model 
which can describe well the elementary electroweak pion production on a single nucleon in the considered energy 
region. Such a model has been developed by Sato and Lee|0] (SL) and is used in this work. The SL model also 
provides accurate pion-nucleon scattering amplitudes which are derived in a way consistent with the pion production 
amplitudes from the same Lagrangian. The 7zN scattering amplitudes are important ingredients in our calculations of 
an optical potential which describes the essential final pion-nucleus interactions. This optical potential takes account 
of the most important medium modifications of the A-propagation within the well-established A-hole model of pion- 
nucleus scattering. Therefore, a combined use of the SL model and the A-hole model, which seems promising for 
calculating pion production in nuclei, provide us with a consistent microscopic description of the pion scattering 
and the electroweak pion production mechanisms in nuclei. This consistency is an appealing point of our approach 
over the previous microscopic models for the coherent pion productions, and enables us to: (i) fix all free parameters 
in our model using the pion-nucleus scattering data; (ii) perform parameter-free calculations for the coherent pion 
productions; (iii) test the reliability of our approach using data for the photo-production process. Another important 
point to be noted is that we take care of the non-local effects on the in-medium A-propagation whose possibly large 
effect was pointed out recently 1 7] . For the neutrino-induced coherent pion production, no previous calculations (either 
in the PCAC-based or microscopic model) have included this effect. 

Our calculations proceed as follows. We first employ the on- and off-shell nN scattering amplitudes generated from 
the SL model to construct a pion-nucleus optical potential within the A-hoIe modeliH. The parameters of the model are 



determined by fitting the pion-nucleus scattering data. By using the pion-nucleus scattering wave functions generated 
from the constructed optical potential and the electroweak pion production amplitudes generated from the same SL 
model, we can calculate the coherent pion production cross sections. Our second step is to establish the reliability of 
our model by comparing the predicted pion photo-production cross sections with the available data. We then proceed 
to calculate the neutrino-induced coherent pion production cross sections. In the following sections, we explain our 
calculational framework, show some selected numerical results, and then conclude. 

A DYNAMICAL MODEL 

We briefly explain how we calculate the amplitude Ai,^j^, for the coherent pion production off a nucleus in its ground 
state (f) induced by an external current denoted by A. Our starting point is a set of elementary amplitudes generated 
from the SL model, a^^^j^^ = aa + a,,,-, where aa (««, ) denotes the resonant (non-resonant) part. The resonant part 
can be written as oa = N/D{W), where D{W) ~W — m^ — '^{W) with W, ;«^, E denoting the nN invariant mass, the 
bare A mass, and the A self energy, respectively. Within the A-hole modeltSJ, we can include the medium effects by 
modifying the A propagator and obtain the following expression of nuclear amplitude Ai,^j^, 

where q and k (q and k) are the momenta for the external current and the pion in the pion-nucleus (pion-nucleon) 
center-of-mass frame; — + q — p^ + ^- A single nucleon wave function in the initial [final] nucleus is denoted 
by WjiPN) [VO(A'a')] ^i'^h the index j specifying a nucleon orbit including the isospin state. The summation (^y) 
is taken over the occupied states of the nucleus. The medium effects on the A-propagator are described by the A 
Hamiltonian {H\), the Pauli-correction to the A self-energy (Zpauii), and the so-called spreading potential (2spr)- The 
spreading potential has the central and spin-orbit parts with their strengths determined by two free complex parameters 
in our model. The total energy of the pion-nucleus system is E +AmN (A: mass number). The integration over in 
Eq.(l) can be done analytically by fixing the p^-dependence of the function in the square bracket to an effective 
constant value, except for the kinetic term in H/^ which is is the source of the non-locality of the A propagation in 
nuclei. We then can sum over the nucleon states to obtain an expression which depends on the nuclear density; the 
nuclear density has been determined experimentally. 

A full transition amplitude is obtained by convoluting the amplitude in Eq. ^ with the pion scattering waves. 
In this work, the pion scattering wave functions are obtained with an optical potential {U:,it). Our starting point of 
constructing Uj^t is the kN scattering amplitude generated from the SL model T^^ = fA + tnn where fA and f„,- are the 
resonant and non-resonant amplitudes, respectively. Consistently with Eq. ([1]), the A-propagator of the resonant part 
fA — FjcnaFjcna/ D{W) is modified by the same procedure. The resulting optical potential is given by 

The last term of the above expression is a phenomenological term proportional to the square of the nuclear density. 
It describes the absorption of s- and p-wave pions by a pair of nucleons and hence is determined by two complex 
couplings Cj and Cp. Thus we have four free parameters in our model: Cj and Cp, and two parameters specifying 
the strength of the spreading potential which are needed to calculate Espr in Eqs.(l)-(2). We determine these four 
parameters by fitting to the pion-nucleus scattering data. We then can perform parameter-free calculations of the 
coherent electroweak pion production cross sections. 

RESULTS 

Pion-nucleus scattering 

We determine the four free parameters of our model by fitting to the available pion-nucleus scattering data in the 
50 MeV <Tji< 300 MeV region. (Tj^ is the pion kinetic energy in the laboratory frame.) Two sample results from 
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FIGURE 1. (a) Total cross sections for — '^C scattering. Tlie solid curve is obtained with our full calculation, while the dashed 
curve is obtained without the spreading potential, (b) Tt^ — '^C elastic differential cross sections at Tjz = 180 MeV. The solid curve 
is obtained with our full calculation. In both figures, the data are fromRef. d. 

our fits to the n — '^C scattering data are siiown in Fig. [T] In the left-hand side of Fig. [H we see that the total cross 
sections for tt^-'^C scattering as a function of Tji can be reproduced very well by our full calculations (solid curve). 
For a comparison, the results obtained without including the spreading potential are shown in the dashed curve. We 
observe a large reduction from the dashed to solid curves, indicating the importance of the strong pion absorption 
simulated by the spreading potential. In the right hand side of Fig. [T] we see that our model can also describe very well 
the differential cross sections. Overall, our model can satisfactorily reproduce the data for both the total and elastic 
cross sections in the considered 50 MeV < T^r < 300 MeV energy region. 

With the four parameters of our model determined in the fits to the pion-nucleus scattering, we are now in a position 
to perform parameter-free calculations of the coherent electroweak pion production cross sections. 

Coherent pion photo-production 

The photo-production process, for which extensive data are available, provides a good testing ground for our 
approach. In the left-hand-side of Fig. |2]we compare the existing data with our results for the differential cross sections 
of 7+ ^ tt" + i^c^^,. . The solid curve is the result of our full calculation. The long-dashed curve is obtained 

without including the pion-nucleus final state interaction (FSI) and without including the the medium effects (Epauii 
and Lspr) on the A-propagation of Eq.(l). When the the medium effects on the A propagation in Eq.(l) is included, 
we obtain the short-dash curve. By comparing these three curves, it is clear that the medium effects are quite sizable, 
and they play an important role in bringing the calculated differential cross sections in agreement with the data. The 
good general agreement seen in Fig.|2|a) indicates the basic soundness of the method we have used in determining the 
spreading potential. It is true that, in the relatively larger angle region, there are noticeable discrepancies between the 
results of our full calculation and the data. We remark however that, as noted in Ref. ifToll . the data in this region are 
likely to be substantially contaminated by incoherent processes in which the final nucleus is in its low-lying excited 
states. 

As a comparison, the dash-dotted curve in left hand side of Figure |2] is obtained when the non-resonant pion 
production amplitud43 is set to zero. Clearly, the non-resonant production is substantial and can not be neglected. 
In general, we observe that, even near the resonance energy, the contribution from the non-resonant mechanism is 
quite significant. This is partly because the resonant contribution is considerably suppressed by pion absorption (the 
spreading potential) and the non-local effect of A propagation (A kinetic term). 

Through the comparison of our numerical results with the pion photo-production data, we have established the 



' In the SL model, the resonant amplitude itself contains the non-resonant mechanisms. We refer to purely non-resonant amplitudes as "non- 
resonant amplitudes". 
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FIGURE 2. (a) Differential cross sections for y- 



■ 1^^ + '^Ci; .5.. The solid line represents the result of the full calculation. 



The long-dashed line is obtained without the FSI and without the medium effects on the A-propagation, while the short-dash line is 
obtained with the medium effects on the A included. The dash-dotted line corresponds to a case in which the pion production 
operator includes only the A mechanism. The data are from Ref. tlQI . (b) The £v -dependence of the total cross section for 
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FIGURE 3. (a) The pion momentum distribution for Vp + C^,.s. /x^+;r++ Cg..s, at iiy = 1 GeV; pji is the pion momentum 
in the laboratory frame. The use of the solid, dashed, dotted and dash-dotted lines follows the same convention as in Fig.|2ta). (b) 
The flux-convoluted rj -distribution for v + '^Cg.^, ^ v + + '^Cg.j, obtained in our full calculation. The neutrino flux is taken 
from MiniBooNE fTTIl . Also shown is the Monte Carlo result from MiiriBooNE Q] rescaled. 



reliability of our approach based on a combined use of the SL model and the A-hole model. Thus we proceed to apply 
the same approach to make predictions for cross sections of the neutrino-induced coherent pion production on '^C. 



Neutrino -induced coherent pion production 

We now present our predictions of the neutrino-induced coherent pion production on the '^C target. We consider 
the CC and NC processes induced by a neutrino or an anti-neutrino. In the right-hand-side of Figure |2l we show the 
predicted total cross sections for these processes as functions of the incident neutrino (anti-neutrino) energy Ey in 
the laboratory system. It is seen that, for higher incident energies, the ratio Occ/<^nc approaches 2, a value expected 
from the isospin factor. For lower incident energies {Ey < 500 MeV), however, a^c is larger than acc, reflecting 
the fact that the phase space for the CC process is reduced significantly by the muon mass. To compare our results 



with data, we need to evaluate the total cross sections averaged over the neutrino fluxes that pertain to the relevant 
experiments. We choose to use the fluxes up to iiv < 2 GeV and neglect the fluxes beyond that limit. We also need 
to consider the experimental setting such as kinematical cuts. A K2K experiment |1] reports the upper limit for the 
neutrino CC coherent pion production, C7k2k < 7.7 x lO^'^^'cm^. This upper limit corresponds to events satisfying the 
muon momentum cut, > 450 MeV and the cut on the square of the momentum transfer, 2^^, < 0.1 GeV^; Q^^ is 
the momentum transfer reconstructed with an assumption of the quasi-free kinematics. Therefore, we also calculate 
the total cross section with these cuts and then convolute with the flux reported by the K2K experiment 1 12]. We obtain 
•'are = 6.3 X lO^^'Vm^ which is consistent with the upper limit from K2K. 

For the NC process, we use the flux reported by MiniBooNE in Ref. fill] and arrive at = 2.8 x 10^'*''cm^; no 
kinematical cut is necessary in this case. This is to be compared with aiviiniBooNE = 7.7 ± 1 .6 ± 3.6 x lO^^'^cm^ given 
in Ref. |13]. Our result is consistent with the empirical value within the large experimental errors, even though the 
theoretical value is rather visibly smaller than the empirical central value. It is to be noted however that Ref. 1. 1 3.1 is a 
preliminary report, and that, as discussed in great detail in Ref. 13, OMiniBooNE may be overestimated due to the use 
of the Rein-Sehgal (RS) modelfis'] in the analysis. 

Next, we present our result for the pion momentum spectrum for the CC process aXEy = \ GeV in Fig. |3la). The 
importance of the medium effects manifests itself here in the same manner as in the photo-process [Fig.|2ja)]. In the 
A region, strong pion absorption is seen to reduce the cross sections significantly, and FSI shifts the peak position. 
The dash-dotted line corresponds to a case in which the pion production operator contains only the A mechanism 
(without non-resonant contributions), while the pion optical potential is kept unchanged. We note that, at /Sy = 1 GeV, 
the dash-dotted line corresponds to 82% of the solid line, and for = 0.5 GeV, 64%. The result indicates that the 
non-resonant components in our model play a significant role in the coherent pion production; their role is particularly 
important for < 0.5 GeV. This characteristic feature of our model should be contrasted with the fact that non- 
resonant mechanisms play essentially no role in any of the previous microscopic calculations for neutrino-induced 
coherent pion production. 

We examine the effects of the non-locality of A-propagation in nuclei arising from the A kinetic term. For this 
purpose, we introduce a quantity ^{Ey) = G{Ey)/ aic,ca\{Ev), where a (diocai) is calculated with (without) the A 
kinetic term. This subject has been studied in Ref. [7] which included the A-mechanism only, without FSI or medium 
modifications of the A. They showed that the cross section is changed by ^ < 0.5, ^ 0.6, < 1 at£'v=0.5, 1, 1.5 GeV. 
When we use the free A-propagator without FSI, which contains the physics similar to that of Ref. [7], we find that the 
non-locality change the cross sections by = 0.41, 0.74, 0.88 at Ey = 0.5, 1, 1.5 GeV. This result is fairly close to 
that obtained in Ref. |7]. Next we use the A-propagator in Eq. ([TJ but without the major medium effects (Espr, Xlpauii, 
and small effects), and then find that the effect of the non-locality is somewhat moderated: ^ = 0.54, 1.03 and 1.14 
at Ey = 0.5, 1.0 and 1.5 GeV. Including the non-resonant amplitudes, medium effects on the A and FSI, the effect of 
the non-locality is ^ = 0.34, 0.73, 0.84 at Ey = 0.5, 1, 1.5 GeV for the full calculation, showing that the non-locality 
is still important. Even though the non-locality might be partly accounted for with the use of the spreading potential 
fitted to data, considering its importance, it is preferable to take it into account explicitly. 

The MiniBooNE analysis of NC data used the tj -distribution [ Tj = Eji{l — cos 9n)] to distinguish coherent pion 
production from the other processes contributing to the TT^-production events. To this end, MiniBooNE used the 
"shape" of the tj -distribution obtained from the RS model [15] with the momentum reweighting function applied. 
However, a microscopic calculation in Ref. [14j was found to give an T] -distribution appreciably different from that 
obtained in the RS model, and the authors of Ref. [[Till pointed out a possibility that the MiniBooNE might have 
substantially overestimated the NC events. Thus it is interesting to compare our result with the T] -distribution from 
MiniBooNE. Figure [3lb) shows the "average" t] -distribution, resulting from the convoluting of the 77 -distribution 
obtained in our present calculation with the MiniBooNE neutrino flux jl 111 . The figure also presents the MiniBooNE 
Monte Carlo results {cf. Fig. 3b of Ref. ]3]), arbitrarily rescaled to match the theoretical curve at 77 = 0.005 GeV. We 
remark that the 77 -distribution we have obtained is fairly close to that given in Ref. iflill . and thus also arrive at the 
conclusion that it is possible that MiniBooNE might have substantially overestimated the NC events. 

DISCUSSION AND SUMMARY 

Recent theoretical efforts on the neutrino-induced coherent pion production are indeed active, and it is interesting 
to see a numerical comparison of those calculations. Such a comparison has been presented at this conferencelfl6|. 
Because of the lack of data, it is impossible to rate the models from the comparison. In the following, however, we 
argue that our present calculation would have advantages over the existing calculations. Our model provides, thanks 



to the consistency built in the SL model, a consistent description of the pion-nucleus scattering and the photo- or 
neutrino-induced coherent pion productions. Because of the consistency, we were able to fix the parameters using the 
data for the pion-nucleus scattering and to predict the coherent processes. It is also the consistency which enables us 
to test our approach using the existing data for the photo-process. No other microscopic calculation for the neutrino- 
induced coherent pion production has been checked to this level. Besides, all the other previous calculations lack the 
non-locality of A-propagation whose importance was examined in the previous section. Further remark is concerned 
with FSI. Our optical potential is made to reproduce the data for the pion-nucleus scattering in the range of 50 MeV < 
Tji ^ 300 MeV, covering the most important A region. On the other hand, in the other microscopic calculations, their 
optical potential may not be as accurate as ours. For example, let us look into the model of Amaro et al. 1 14], which 
is the most sophisticated among the existing microscopic models for neutrino-induced coherent pion production. The 
parameters in their optical potential were fitted to the data for the binding energies of pionic atoms. Because, in the 
present context, the optical potential should work well over relatively wide energy region around the A region, the 
optical potential used in Ref. ||14|] may contain a questionable aspect. Finally we mention that it has become fairly 
clear that the RS model, the prominent PCAC -based model, does not work reasonably for Ey <2 GeV. This point was 
made clear through the detailed study of the relation between the RS model and the microscopic model in Ref. 1 14]. 
Although some improvements on the RS model has been proposed iUtI ITsi fl9ll . the applicability of a PCAC-based 
model to the relevant energy region is currently a subject of controversy. 

In summary, we have developed a microscopic model based on a combination of the SL and the A-hole models. The 
model provides a consistent description of the pion-nucleus scattering and the coherent pion production processes. 
Utilizing the consistency, we can fix the free parameters in our model using the pion-nucleus scattering data, and then 
give parameter-free predictions for the coherent processes. After testing the reliability of our approach using the data 
for the photo-process, we calculate and present results for the neutrino-induced coherent pion production. 
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